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In the solid-state circuit QED system and based on the homodyne measurement in dispersive regime, we 
demonstrate that a homodyne-current-based feedback can create and stabilize highly entangled two-qubit states 
in the presence of moderate noisy environment. Particularly, we present an extended analysis for the current- 
based Markovian feedback, which leads to an improved iiltered-current-based feedback scheme. We show that 
this is essential for us to achieve the desirable control effect in present system. 
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I. INTRODUCTION 



The circuit quantum electrodynamics (QED) lH-yl], a solid- 
state analog of the conventional quantum optics cavity QED, 
is a promising solid-state quantum computing architecture. 
This architecture couples superconducting electronic circuit 
elements, which serve as the qubits, to harmonic oscillator 
modes of a microwave resonator, which serve as a "quantum 
bus" that mediates inter-qubit coupling and facilitates quan- 
tum measurement for the qubit state. At the early stage of ex- 
periments, strong coupling between the qubits and the quan- 
tum bus in the circuit QED system should make it excellent 
platform for quantum control study. 

In particular, quantum measurement in this system can be 
carried out by operating in the dispersive limit (where the de- 
tuning between the resonator and the qubit is much larger than 
their coupling strength). In this limit, the interaction induces 
a qubit state dependent frequency shift on the resonator By 
measuring the resonator output voltage with a homodyne mea- 
surement, information about the qubit state is obtained. In 
our present work, based on this homodyne measurement to- 
gether with the flexibility/advantage that feedback can be ap- 
plied to either the qubit or the microwave resonator mode, we 
will illustrate how highly entangled state can be created and 
stabilized from an initially separable state, by an appropriate 
feedback control. It has been well known that quantum entan- 
glement is one of the key ingredients in order to realize quan- 
tum computation and quantum information processing. Very 
recently, in the circuit QED system, interesting ideas were 
proposed to probabilistically create entangled states by means 
of the homodyne measurement alone |0], following the gen- 
eral idea that measurement can be used as a nondeterministic 
means of preparing quantum states that are otherwise diffi- 
cult to obtain. However, besides the drawback of probabilistic 
nature, this measurement alone approach has no ability to sta- 
bilize/protect the obtained entangled state. Generally speak- 
ing, entanglement degradation through uncontrolled coupling 
with the environment remains a major obstacle in practice |J5|- 
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131, which would limit the lifetime of entangled states and de- 
mand efficient schemes to protect them. In this context we 
remind that, owing to the remarkable progress in theory and 
particularly in experiments on the real-time monitoring and 
manipulation of individual quantum systems |@|, the quantum 
feedback control technique may emerge as a natural possible 
route to develop strategies to prepare entangled states and pre- 
vent their deterioration. 



In Ref. iQl, for the optical cavity QED system, the steady- 
state of two qubits (two-level atoms) interacting simultane- 
ously with a driving laser was shown to be entangled, achiev- 
ing a concurrence of 0. 11 without measurement or feedback. 
A later work, involving the use of a homodyne-mediated di- 
rect feedback modulation of the laser that drives the cavity 
mode llioll . extended this study and showed that one can in- 
crease the steady-state entanglement, i.e., the amount of the 
maximum steady-state concurrence can be increased from 
0.11 to 0.31. In principle, maximally entangled states could 
be achieved l\M. by the use of Bayesian, or state estima- 
tion, feedback lll2ll . This improvement, however, comes at 
the price of increasing the experimental complexity, due to 
the challenging need for a real-time estimation of the quan- 
tum state. More recently, rather than the homodyne-mediated, 
a jump-based direct feedback is demonstrated to allow the ro- 
bust production of highly entangled states in the optical cavity 
system 1 1 3il . For the solid-state circuit QED, however, the sin- 
gle photon counting at microwave frequencies is currently not 
possible, despite a couple of proposals and efforts out there 
for doing it liRITsIl 



In this paper we base the feedback on the more available ho- 
modyne measurement scheme [160. Particularly, in Ref. iflTll . 
detailed analysis for measurement scheme to reach the quan- 
tum limit was presented, and a deterministic creation of entan- 
glement by using the measurement based feedback was illus- 
trated in the absence of environmental decoherence. In present 
work we extend the study in Ref. 111711 by taking into account 
the environmental influence. We will show that all the four 
Bell states can be deterministically created and stabilized in 
the presence of moderate noisy environment, with high qual- 
ity compared to the existing results mentioned above. 
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FIG. 1: Schematic diagram of circuit QED together with a mi- 
crowave transmission measurement and a measurement-current- 
based feedback loop. The Cooper-pair box qubits are fabricated in- 
side a superconducting transmission-line resonator and are capaci- 
tively coupled to the voltage standing wave. 



II. MODEL AND FORMALISM 

In Fig. 1 we illustrate the schematic setup of the solid-state 
superconducting circuit-QED system, together with the mea- 
surement and feedback control idea. The central section of su- 
perconducting coplanar waveguide plays the role of the cavity 
and the superconducting qubits play the roles of the atoms. 
The superconducting qubits are coupled to a one-dimensional 
transmission line (IDTL) cavity which acts as a simple har- 
monic oscillator. The qubits, the IDTL cavity and their mu- 
tual coupling can be well described by the Jaynes-Cummings 
Hamiltonian: 



H - a>ra^a + &(a ' +a)+ y 



Q, 



—o-~. + gj{o-ja^ +a-ja) 



(1) 



Here, the operators crT(cr|) and fl(a') are, respectively, the 
lowering (raising) operators for the yth qubit and the cavity 
photons. cOr is the frequency of the cavity photon, and Q; and 
gj are the jth qubit transition energy and coupling strength to 
the cavity photon. For simplicity and for the purpose to be 
clear later, in this work we assume that Qi = Q2 = fl, and 
g\ - -gi - 8- This implies that we assume two identical 
qubits located in the cavity at places with the maximum field 
amplitude and with a half-wavelength separation, as schemat- 
ically shown in Fig. 1. The £-term in Eq. ([T]) stands for a 
microwave driving to the cavity that is employed here for the 
task of measurement. More explicitly, £ = ee"'""''+c.c, where 
the frequency can differ from the cavity photon frequency, i.e., 
Ar = COr - Wm ^ 0. In conccm with the qubit-cavity coupling, 
we focus on the dispersive regime JlHst], which corresponds 
to an energy detuning, A = w, - Q, much larger than g. In 
this limit, the canonical transformation, Heff - U^ HU, where 
U - exp[2j /lj(acr| - a^cr^)] with Aj = gj/^, yields (in the 
rotating frame with the driving frequency (d„,) 






(2) 



<j\ + o\. 



Here, we defined ;t' = g^ 1^ and J, 

In the circuit-QED system, the measurement is imple- 
mented via a homodyne-type detection of the transmitted mi- 
crowave photons as schematically shown in Fig. 1. For the 



composite system of the qubits plus the IDTL cavity, the leak- 
age of photons is described by a Lindblad term KD[a]p in 
master equation, where k is the leakage rate and the Lind- 
blad superoperator acting on the reduced density matrix p 
is defined by D[a]p = apa' - \{a'a,p]. However, con- 
ditioned on the output (homodyne) current, i.e., /hom(0 = 
K{a + a^)c(t) + \[K^(f), there will be an additional unravelling 
term in the conditional master equation, 'J-([a]pc^{t). Here, 
((■ ■ ■ ))c(t) = Tr[(- ■ ■ )pcit)] with pc(0 the conditional density 
matrix, and 'J-([a]pc = ap^ + PcCi^ — Tr[(a + a^)pc\pc- And, 
the quantum-jump related stochastic nature is characterized 
by ^(f), a Gaussian white noise with properties of £[^(0] = 
and E{^(t)^{t')\ - 6(t — t'), where E[- ■ ■] means an ensem- 
ble average over realizations of the noise. Furthermore, in 
this work, we are interested in the regime of strongly damped 
cavity, which enables us to adiabatically eliminate the cavity 
degree of freedom flSll . Qualitatively, from an observation on 
the effective coupling x<^^<^Jz in Eq. (|2]i, we can expect: the 
measurement-backaction induced dephasing term ~ D[Jz]pc, 
the unravelling term ~ '7Y[7;]pc^(0, and the homodyne current 
homit) ~ (/-)e(0 +^(t)- Indeed, following the standard proce- 
dures of adiabatic elimination jTSlI , one obtains the quantum 
trajectory equation (QTE) involving only the qubit degree of 
freedom, in a rotating frame with respect to the qubit Hamil- 
tonian which reads 

p, - -ix\af[J„pc] + Yj rPkJK + X ^^['^•H 



j=h2 



>=1.2 



+ypD[cr- - cT,]p, + ^D[J,]pc + ^•H[7,K^(r). (3) 

In deriving this result, we have used the assumption A,. = 
and Ai = -A2 = A. In Eq. (|3]l, yj and -y^j are the relaxation 
and dephasing rates caused by the surrounding environment, 
while -/p = kA~ is a rate for the collective decay due to the Pur- 
cell effect. Explicitly, the measurement-backaction induced 
dephasing rate F,/ = 8|Qf|y //c, with a = -Hs/k. Finally, 
the information gain rate F,,, in Eq. ([3]) is in general related to 
the backaction dephasing rate through the quantum efficiency, 

77 = F,„/(2rd). 



III. HOMODYNE-MEDIATED FEEDBACK 

A. Feedback Equation in Markovian Limit 

Based on the QTE, one can infer the conditional state pdt). 
In principle, one can then perform a state-estimate based feed- 
back, by an appropriate design for the feedback Hamiltonian. 
This is the so-called Bayesian feedback scheme. Another, 
much simpler, scheme is the direct feedback which depends 
linearly on the detection signal, e.g., the homodyne current in 
the present case. Indeed, the former scheme can result in an 
improvement over the direct feedback, yet it comes at the cost 
of an increasing experimental complexity due to the challeng- 
ing need for a real-time estimation of the quantum state. In 
what follows, we focus on the direct current-based feedback 
scheme. In general, we denote the feedback Hamiltonian as 



//fb(0 - Ihomit - t)/^, where t stands for a possible time delay 
of the feedback, while the homodyne current can be expressed 
as Ihomit) - ^fr^{Jz)c(t) + f(0 after the adiabatic elimination 
of the cavity degree of freedom. Performing this feedback 
leads to state change according to [p^ (f)]fb - homit " TYKpdf) 
, where 'Kpdt) = —i[F,pc{f)\. In the Markovian limit by as- 
suming T = 0, an appropriate combination of this equation 
with the QTE © yields 

p, = -ix\a\^[J,,pc] + ^D[J,]pc 

+ 2 jjD[crJ]p, + Yj ^^D[crj]p, 



J=l,2 



j=l.2 



Vrl, 



+ypD[o-i - a-^jpc + —r-'KiJrPc + PcJd 



+ \'K^Pc + [^'Hm+'K\pm- 



(4) 



Originally, this type of equation was obtained by Wiseman 
and Milburn 111911 . by a careful interpretation to the feedback 
as an Ito- or Stratonovich-type stochastic action. Interest- 
ingly, we mention here an alternate (and equivalent) deriva- 
tion for this result. Since even in the Markovian limit the 
feedback can be applied only after the measurement out- 
come, in terms of discretized time interval we thus have 
Pc{t + dt) - e''^''*pc(f + dt), where pdt + dt) is the state af- 
ter inferring from the measurement record, and the differen- 
tial superoperator reads d'K{t) - dlhomitYK. Here, dibomit) - 
y/T7,{J-)c(t)dt + dW(t), while dW(t) = ^(t)dt is the Wiener 
increment that has the statistical properties of E[dW(t)] = 
and E[dW{t)dW(s)] = 6{t- s)dt. Expanding this instantly fast 
feedback to second order and keeping the infinitesimal incre- 
ment to 0{df) (noting that [(iW(0]^ - dt), one then obtains 
Eq.©. 



Entanglement Creation and Stabilization: Preliminary 
Result 



From now on we specify our study by using feedback to 
create two-qubit entanglement, first based on Eq. (HI, then on 
an improved scheme. In particular, we will focus on feed- 
back protection against the spontaneous emission of qubits, 
while assuming the dephasing rates y^j negligibly small for 
the transmon-type qubit ^. For two qubits, there are four 
maximally entangled states, i.e., the Bell states; \^'±} - 
^(|00>+|ll»,and|(D±> = ^(|01)+|10)). In what follows we 
only detail our analysis for IO+), and remain the others as brief 
discussion in the final concluding section. Initially, we assume 
a separable state for the two qubits, (|0) + |l»i ® (|0> + |1»2 = 
(|00) + 111)) + (|10) + |01)), which can be prepared easily by 
separate single bit rotations. Within the scheme after adi- 
abatic elimination of the cavity degree of freedom, we see 
that the current {Jz)c is zero for qubit state |10) + |01), and 
nonzero otherwise. Then, we expect that the current-based 
feedback would force the state towards the entangled state 
|10) + |01), if we exert a 7j-type flipping, i.e., apply a feed- 



CO.5 




C0.5 



0.0 



f\ 






■1 


I 1 Y = o.oir^ 


- 


r 

IF i 1 i 1 1 1 i 1 1 


, (c)' 



10 15 20 



25 



30 



FIG. 2; (a): Average concurrence over 500 trajectories based on 
Eq. ^ and a direct current feedback //fb = j(/hom(0-/j:- (b) and (c): 
Concurrence obtained by a feedback Hfb = u{J-)c(t)Jx, i.e., remov- 
ing the noise in the homodyne current. In (b) the average and in (c) 
the individual quantum trajectory concurrences are illustrated, both 
showing considerable improvement over the result in (a). Parame- 
ters: 7,, = Fj, M = 0. 1 in (a) and u = 1.0 in (b) and (c). Here, the 
choice of different feedback strength u = O.l in (a) is from a rough 
numerical optimization, since an increase of u over this value will 
give poorer result instead. 



back Hfb{t) - uIbom{t)Jx, where u characterizes the feedback 
strength over the current-based modulation. 

In Fig. 2(a), we illustrate the control result based on Eq. (01, 
where the average concurrence over 500 trajectories is shown. 
Note that the concurrence, which is particularly useful to char- 
acterize two-qubit entanglement in a mixed state[2(|l, is an 
effective measure in present context for the creation and sta- 
bilization of !<!'+), which is in good consistence with the state 
fidelity, F{t) = Tr[|(D+)<0+|p(f)]. Compared to Refs. lH and 
lUOll . where the achieved concurrences are, respectively, 0.11 
and 0.31, we see that the result in Fig. 2(a) is good, but not 
optimal for the following reasons. First, we note that this 
improvement over the result in Ref. lllOll is largely owing to 
the different measurement scheme here, which in the disper- 
sive regime is in fact a continuous quantum non-demolition 
(QND) measurement. That is, here the measurement is a J.- 
type, while in Ref. jlOll it is a yj-type. For current-based feed- 
back, as heuristically discussed above, the homodyne current 
dihomit) contains, not only the useful signal {J:}c{t), but also 
the harmful noise (i.e. the dW{t) term). To this point, it is im- 
portant to understand the distinct role of the dW{t) term in the 
homodyne current, when using it to update the quantum state 
versus to perform the feedback. In doing the former, it is infor- 
mative; while in doing the latter, it is useless and harmful. In 
the current-based Markovian feedback equation, i.e., Eq. (|4|i, 
the two Wiener increments were equated and were mixed up 
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FIG. 3: (a) Average concurrence over 500 trajectories, and (b) the 
concurrence of three representative individual trajectories, by apply- 
ing a filtered-current-based feedback control. Parameters: y,, = Fj, 
u = 10, 7f, = 0.006F,,, and T = 2000 * dt. 



in the final form ||19ll . We thus understand that it may lead 
to inefficient result sometimes. As a quantitative support to 
this reasoning, in Fig. 2(b) and (c), we show the result by 
adding a feedback, Hfb - u{J~)c{t)Jx, straightforwardly into 
the measurement record conditioned quantum trajectory equa- 
tion. We remind here that the noise in the homodyne current 
was removed. In Fig. 2(b), the average concurrence by aver- 
aging over 500 quantum trajectories is plotted, while in Fig. 
2 (c) the concurrence of individual quantum trajectories is il- 
lustrated. Noticeably, essential improvement over the result in 
Fig. 2(a) is achieved. 



C. FUtered-Current-Based Feedback: Improved Result 

As a matter of fact, the above "improved" feedback, //ft, - 
u{Jz)c{t)Jy:, is a State-estimation feedback, since the {J~)c{t) is 
known only after knowing the state pdt). However, guided 
by it, we can refine the noisy homodyne cuiTent. Following 
Ref. iflTll . we first low -pass filter the measurement signal over 
a small time window, R{f) — jj J e"^"*'"'"'(//hom('''), where 
the factor N normalizes the maximum of the smoothed sig- 
nal R{t) to unity. Through this filtering procedure, we get 
in fact a crude but efficient estimate of (7,)f.(f). Then, we 
condition the feedback with a power of the filtered measure- 
ment signal, which would reduce the noise further in the es- 
timate. That is, we finally design a feedback Hamiltonian, 
Hfhit) - uR{t)^J-(, adding directly into the measurement con- 
ditioned quantum trajectory equation, where P is the power to 
which the smoothed signal is raised. In Fig. 3 we display the 
(a) average and (b) individual trajectory concurrences, based 
on this filtered-current feedback control. We find that the re- 
sult is satisfactory, being comparable to that in Fig. 2(b) and 
(c). 



FIG. 4: Environment caused deterioration of the Bell state |cD+> in 
the absence of feedback protection. In (a) the average fidelity (over 
1000 trajectories) and in (b) the fidelity of a representative trajectory 
are presented (blue curves), while the result with feedback is plotted 
here for comparison (red curves). The feedback parameters are the 
same as in Fig. 3. 



IV. ENTANGLEMENT DETERIORATION IN THE 
ABSENCE OF FEEDBACK 



In addition to the deterministic creation of entanglement 
achieved here, we remark that the other advantage of feed- 
back is its ability to stabilize the entanglement. We illustrate 
this point more explicitly by the following. Given the entan- 
gled state |<1>+) having been achieved, for instance, probabilis- 
tically by the measurement alone approach as proposed in Ref. 
H, in Fig. 4 we plot the later fate of this state in terms of 
its fidelity F(t) - Tr[(|O+>(<I)+|)p(0], in the absence of pro- 
tection with feedback. In Fig. 4(a) we show the average fi- 
delity for the ensemble averaged state, which is found to de- 
cay with time under the influence of environment. By altering 
the coupling strength to environment, this decay may turn to 
the remarkable phenomena of sudden death of entanglement 
151], which were discussed under a broad class of dissipation 
models. Meanwhile, in Fig. 4(b), the fidelity of representative 
individual trajectory is shown. Interestingly, the sudden jump 
of fidelity to zero is observed, which indicates a sudden disap- 
pearance of entanglement in our present case. We stress that 
this is a real "sudden death" of entanglement, in the sense of 
an instantaneous disappearance of entanglement, in contrast 
to the gradual death (i.e., with finite lifetime) of the ensemble 
averaged state ||5[]. Anyhow, we may conclude, that even the 
entangled state (e.g. |0+)) has been probabilistically created 
|4[|, it will completely die, probably after some short surviving 
time, but almost certainly with the increase of time, in the ab- 
sence of feedback protection. Nevertheless, the feedback we 
discussed above can prevent this from happening. 



V. DISCUSSION AND SUMMARY 

So far we have examined a feedback scheme of determin- 
istic creation and stabilization of entanglement. As a specific 
example, we considered the state |<1)+) = -L(|01) + |10)) in 
detail. In practice, the need of "/v (=crf + cr^y'-type feedback 
can be realized by a modulation of the qubit Hamiltonian pa- 
rameters, as schematically indicated in Fig. 1, or more sim- 
ply, by modulating the microwave driving strength. For the 
latter scheme, besides the measurement driving, one can ap- 
ply an additional control microwave driving, with amplitude 
of fc and frequency in resonance with the qubit transition en- 
ergy. It can be shown that this driving will induce an effective 
Hii = /lef/v on the qubits, if we properly set the qubits so 
that gi - g2 - g, or equivalently, Ai - A2 - A. Notice that 
in this case the Purcell term turns out to be ypD[cri + cr2]p, 
instead of ypD[cri - cr2]p as in Eq. ([3]), making the state \<i>+} 
not invariant under its action. However, one can make this 
term relatively small by appropriately increasing the detuning 
between the cavity photon and the qubit. Moreover, since this 
term only transits |<l)+> to |00>, yet 7,100) - |0+>, we thus ex- 
pect that the applied feedback can well suppress its destructive 
effect. 

Using similar method, one can deterministically create and 
stabilize another entangled state |<1)_> = 4=(|01> - |10)). This 

can be realized by a homodyne-mediated /v-type feedback, 
where J^ = crj - cr^. This feedback can still be implemented 
by modulating either the qubit Hamiltonian parameters or the 
microwave driving. Finally, in concern with the other two Bell 
states, I*}*-!-) = -1=(|00) ± |11)), we can create and stabilize 
them by a simple conversion scheme, by noting that they are 
in fact related to |<l>±> = -^(|01> + |10)) by a single bit n- 
pulse cr^ rotation. Note that direct creation and stabilization 
of I^F-t) is also possible, by alternatively setting g\ - g2 — g 
and an opposite detuning Qi -co,- - cd^ - Q.2- Furthermore, 
an extension to protect arbitrary entangled states, say, |3)-t) - 
c|01) + <i|10) and \^i) = c|00> + <i|ll), is straightforward, by 
simply replacing the feedback operator 7, by an appropriate 
combination of a-\ and crl. 

Finally, we mention that the main problem with doing feed- 
back in circuit QED is the lack of efficient homodyne detec- 
tion. Currently, the way to perform homodyne and heterodyne 
detection is to first amplify the signal before mixing it on a 
nonlinear circuit element of some kind. As a consequence, the 
extra noise added by the amplifier will reduce the quantum ef- 
ficiency and prohibit quantum limited feedback. It seems that 
this situation is to be changed quickly, for instance, by devel- 
oping Josephson parametric amplifiers which can be realized 
in superconducting circuits lllql . In our present theoretical 
study, we did not include the non-unit quantum efficiency into 



the homodyne detection of the field. This treatment is largely 
owing to the following consideration. After adiabatic elimi- 
nation of the cavity photon degree of freedom, the non-unit 
quantum efficiency of homodyne detection will reduce the ef- 
fective information-gain rate F,,, in Eq. (O. This implies an 
emergence of an extra non-unravelling dephasing term in the 
quantum trajectory equation. However, in the context of cre- 
ation and stabilization of |'I'+), this term only results in de- 
phasing among states |00), |11), and |01) + |10). That is, it 
does not destroy the target state IO+). From the feedback prin- 
ciple discussed in Sec. Ill (B), we can also imagine that this 
dephasing term only affects the detailed control dynamics, but 
does not prevent the state towards the target state. Therefore, 
for certain acceptable (not too low) quantum efficiency of the 
homodyne detection, the present scheme of feedback should 
be implementable in circuit QED. This is seemingly a rare ex- 
ample that the control effect does not sensitively depend on 
the measurement efficiency. We have numerically examined 
the above reasoning, for instance, by lowering the quantum 
efficiency to // = 0.8, and found very small change of the en- 
tanglement concurrence. 

To summarize, in the solid-state circuit QED system and 
based on the homodyne measurement in dispersive regime, we 
analyzed the creation and stabilization of the highly entangled 
Bell states by the use of a filtered-current-based feedback. 
Compared to a few previous studies in similar cavity QED 
systems, we (i) eliminate the drawback of probabilistic nature 
and the inability of stabilization 1^, (ii) improve the control 
effect by enhancing the concurrence to values higher than 0.9, 
in regard to Ref. [lOy where the improved concurrence of 0.3 1 
over the former 0.11 was achieved, and (iii) avoid the exper- 
imental difficulty in present system for the jump-based feed- 
back II13I1 . or the complexity for the state-estimation feedback 
d. 
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